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ABSTRACT 

We report on a comprehensive analysis of the kilohertz (Js 300 Hz) quasi-periodic oscilla- 
tions (kHz QPOs) detected from the neutron star low-mass X-ray binary 4U 0614+09 with 
the Rossi X-ray Timing Explorer (RXTE). With a much larger data set than previously an- 
alyzed (all archival data from Februar y 1996 up to October 20 07), we first investigate the 
reality of the 1330 Hz QPO reported bv lvan Straaten et al.l d2000t) . This QPO would be of par- 
ticular interest since it has the highest frequency reported for any source. A thorough analysis 
of the same observation fails to confirm the detection. On the other hand, over our extended 
data set, the highest QPO frequency we measure for the upper kHz QPO is at ~ 1224 Hz; a 
value which is fully consistent with the maximum values observed in similar systems. Second, 
we demonstrate that the frequency dependence of the quality factor (Q = v/Ais) and ampli- 
tude of the lower a nd upper kHz QPOs follow the systematic trends seen in similar systems 
dBarret et alj|2006l) . In particular, 4U 0614+09 shows a drop of the quality factor of the lower 
kHz QPO above ~ 700 Hz. If this is due to an approach to the innermost stable circular orbit, 
it implies a neutron star mass of ~ 1.9 M Q . Finally, when analyzing the data over fixed dura- 
tions, we have found a gap in the frequency distribution of the upper QPO, associated with a 
local minimum of its amplitude. A similar gap is not present in the distribution of the lower 
QPO frequencies, suggesting some cautions when interpreting frequency ratio distributions, 
based on the occurrence of the lower QPO only. 

tivity (see however llvlendezl d2006h for a different interpretation). 
With the availability of more data in the RXTE archive, it is now 
possible to search for a quality factor drop similar to that seen in 
other sources. 



1 INTRODUCTION 



Kilohertz quasi-periodic brightness oscillat ions (kHz QPOs) 
have been report e d from 4U 06 14+09 by iFord et alj <fl997h; 
Ivan Straaten et alj | |200C|. l2002h : iBarret et alj J2006h : iMendej 
d2006h . The peculiar properties of its kHz QPOs motivate the 
present work for three main reasons. First, among kHz QPO 
sources, 4U 0614+09 holds the record for the highest claimed 



QPO frequency, 1330 Hz ( van Straaten et all I20001) . whereas in 
most sources the maximum frequency for the upper kHz QPO lies 
around 1200 Hz. This is of particular importance because it sets the 
most stringent constraints on the mass and radius of t he NS, under 
the a ssumption that 1330 Hz is an orbital frequency ( Miller et al. 
1 1998h . Unfortunately, 4U 0614+09 tends to have a low count rate 
and broad QPOs compared to similar sources, hence its QPOs are 
challenging to characterize. 

Second, IBarret et"al] d2006h have performed a systematic 
study of the quality factor of the lower and upper kHz QPOs in six 
systems: 4U 1636-536, 4U 1608-522, 4U 1735-44, 4U 1728-34, 
4U 1820-303 and 4U 0614+09. Using data available in the RXTE 
archive at the end of 2004, they found that all the sources except 
4U 06 14+09 showed evidence of a drop in the quality factor of their 
lower kHz QPOs at high frequency. For 4U 0614+09 only the rising 
part of the quality factor versus frequency curve was reported. The 
sudden drop is consistent with what is expected if it is produced by 
the approach of an active oscillating region to the innermost stable 
circular orbit (ISCO), a key feature of strong-gravity general rela- 



Third, 4U 0614+09 has recently drawn further attention, af- 
ter the detection of burst osc illations at 414 Hz with the Burst Alert 
Telescope on board SWIFT I Str ohmaver et al.l2008h . The latter fre- 
quency which is likely to be the neutron star spin frequency (v s ) is 
to be compared with the frequency diff erence At/ ~ 320 Hz be - 
tween twin kHz QPOs reported so far ( Ivan Straaten et al]|2000h . 
Clearly the ratio Av/u 3 was not consistent with either 0.5 or 1, 
observed in similar syste ms. This particul ar re sult was consistent 
with re cent suggestions bv lYin et"aL | d2007h and lMendez & Bellonil 
d2007l) that the kHz QPO frequency difference may not have a 
strong connection to the neutron st ar spin frequency in some 
sources (but see also the discussion in lBarret et alJl20o"8T) . A closer 
inspection of the frequency difference of the twin QPOs over a 
much larger data set is therefore needed in light of this result, but 
als o because a significant s catter is present in the values reported 
bv lvan Straaten et alj d2OO0h . 

Here we perform an analysis of more than 1 1 years of RXTE 
data on 4U 0614+09. In § 2 we describe our analysis scheme and 
present our main results. We discuss the implications of our results 
in §3. 
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2 OBSERVATIONS AND DATA ANALYSIS 

We have retrieved from the HEASARC archive science event mode 
data recorded by the RXTE Proportional Counter Array (PCA). 
The data set spans over eleven years from February 26th, 1996 
to October 17th, 2007. We consider segments of continuous ob- 
servation (ObsIDs): 763 ObsIDs were analyzed with a typical 
duration of 3000 seconds. For each ObsID, we have computed 
an average Power Density Spectrum (PDS) with a 1 Hz resolu- 
tion, using events record ed between 2 and 4 keV. The PDS are 
normalized according to ILeahv et al. I <1983h . so that the Poisson 
noise level is constant around 2. The PDS is then blindly searched 
for excess power between 30 Hz and 1400 H z using a scan- 
ning technique, as presented in lBoirin et alj J2000l) . The frequency 
range searched includes the highest QPO frequency reported so far 
dvan Straaten et al]|2000l) . We have also verified that no significant 
excesses were detected between 1400 and 2048 Hz. This justifies 
the use of the 1400 — 2048 Hz range to estimate accurately the 
Poisson noise level in each observation. Each excess (at most the 
2 strongest) is then fitted with a Lorentzian with three free param- 
eters; frequency, full width at half maximum (constrained to range 
from 2 to 1000 Hz), and amplitude (equal to the integrated power of 
the Lorentzian). The Poisson noise level is fitted separately above 
1400 Hz and then frozen when fitting the QPOs. Errors on each pa- 
rameter are computed with A^ 2 = 1. As in previous papers in this 
field, our threshold for QPOs is related to the ratio (hereafter R) 
of the Lorentzian amplitude to its la enxwQ (R was often quoted 
and used as a significance). In this paper, our threshold for R is 
3, meaning that we consider only QPOs for which we can mea- 
sure the power of the Lorentzian with an accuracy of 3<j or more. 
Such a threshold corresponds to a ~ 6a excess power in the PDS 
for a single trial, equivalent to ~ 4a if we account for the num- 
ber of trials of the scanning procedure dvan der ais|[l989T ). Fur- 
thermore, as expected we have found that R positively correlates 
with the single trial significance of the excess power in the PDS. 
It is worth mentioning that the proportionality coefficient is close 
to 2; i.e. the QPO with the highest R ratio ~ 10 corresponds to a 
~ 20cr (single trial) excess power in the associated PDS. The inte- 
grated power of the Lorentzian is then converted into a root mean 
square (RMS), expressed as a fraction of the total source count rate. 
As said above the QPOs in 4U 0614+09 can be broad. In order to 
recover properly the parameters of the QPOs in the low frequency 
end (~ 400 — 700 Hz), th e continuum underneath t he QPO must be 
accounted for. Following lvan Straaten et al. I i2002h . the continuum 
is adjusted with one or two zero centered Lorentzians. An example 
of such a fit is shown in Fig.Q] showing the decomposition of the 
PDS into several components. 

As a result of this systematic analysis, in 210 ObsIDs we 
detected a single QPO, and in 24 ObsIDs we found two simul- 
taneous QPOs (in one ObsID, two QPOs were detected but not 
simultaneously). Two tables summarizing all the QPO detections 
are available in the electronic version of the paper. The remain- 



1 The Lorentzian function used in the fit is Lor(^) = A X W/ (27r) / \[y — 
vq) 2 + (W/2) 2 ], where A is the integrated power of the Lorentzian from 
to oo, W its width and wo its centroid frequency. The fitted function 
is linear in A, an d therefore its error can be computed using A\ 2 (e.g. 
IPress et al.Hl992l) ). The RMS amplitude is a de rived quantity, comp uted as 
RMS=-y/A/S, where S is the source count rate Ivan der Klisll989l) . In this 
paper, we have defined R = A/ 5 A, from which the error on the RMS is 
estimated as 5RMS = 1/2 X RMS X R" 1 after neglecting the term SS/S 
in the derivative of the RMS equation. 
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Figure 1. Illustration of the multi-lorentzian fit of the PDS for the ObsID 
10095-01-01-00 recorded on April 22nd, 1996 at 22:32 pm. The broad zero- 
centered Lorentzian peaking around 200 Hz must be included in the fit to 
recover pr operly the parameters of t he kHz QPO at 600 Hz. This was also 
noticed bv lvan Straaten et"al] j2002l) . who adjusted the low frequency con- 
tinuum with one or two zero centered Lorentzians. 



ing ObsIDs, which contained no QPOs, were removed from the 
subsequent analysis. The quality factor and RMS amplitude versus 
frequency of all detected QPOs are s hown in Fig [ 2] The patterns 
identified in other similar systems by iBarret et"al] 1 20061) are also 
seen for 4U 0614+09. In particular, in the quality factor versus fre- 
quency plot, lower and upper QPOs occupy two distinct regions: 
the quality factor of lower QPOs is larger than the quality factor of 
upper QPOs and the quality factor of the upper QPOs rises steadily 
with frequency (note that the scatter in the quality factor of the 
lower QPOs is in part because it is not corrected for the frequency 
drift within the ObsID). Similarly, there is a clear trend (albeit with 
some scatter) for the RMS amplitude of the upper QPOs to de- 
crease with increasing frequency, whereas the RMS amplitude of 
the lower QPOs increases first to reach a plateau (see below for a 
full description of the amplitude of the lower QPO). Note also that 
the RMS amplitude of the upper QPO, when its frequency is around 
1 150-1250 Hz lies above the extrapolation of the blue dashed line, 
passing through the RMS values measured at lower frequencies. 
Finally, it is worth mentioning that there is a gap in the frequency 
distribution of the upper QPO, between 1020 and 1130 Hz (this 
point will be discussed in more detail below). 



2.1 On the 1330 Hz QPO 

As can be seen from Fig[2j the highest QPO frequency detected in 
our systematic analysis is at about 1220 Hz; hence we do not de- 
tect a ny QPOs at frequencies s imilar to the one reported at 1330 
Hz by (van Straaten et all d200Lt) . For the 1330 Hz QPO, note that 
our detection threshold (R — 3) corresponds to a single trial sig- 
nificance of 6a, much larger than the 3.5<r single trial significance 
reported by van Straaten et a Still, we have repeated the 

analysis of Ivan Straaten et al.1 1 |2000|) for the ObsID 40030-01-04- 
00 in which the later QPO was reported (considering events from 
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Figure 2. Quality factor (left) and RMS amplitude (right) versus frequency of all QPOs detected in the 2^10 keV range. Each point represents the average 
over one ObsID. Red and blue filled circles are respectively for lower and upper twin QPOs. Black filled squares with green error bars are for single detected 
QPOs, identified as lower QPOs. Black filled circles with green error bars are for single detected QPOs identified as upper QPOs. 



5 keV to 97 keV). The strongest excess we could fit was at v = 
1328.4±26.5Hz,FWHM = 46.2±70.2Hz, RMS = 4.6±1.7%, 
hence with a rati o R of 1.3, a value far to o small to claim a de- 
tection. Note that van Straaten et ail feood) reported a larger ratio 
R ~ 2.8; a value closer to our threshold, but which we failed to 
reproduce (see Figure(3](. 

We have tried to optimize the energy band over which the PDS 
is computed to determine whether the significance of the above ex- 
cess could be increased. By looking at the count spectrum of the 
source, it dominates over the background between 2 keV and 20 
keV. We have thus computed a PDS of this ObsID using only events 
from 2 to 20 keV. The strongest excess of the PDS is no longer 
around 1330 Hz (a R ~ 1 excess exists at 966 Hz). By initializ- 
ing th e parameters of the fit with the values of Ivan Straaten et all 
( feOOOh . the Lorentzian parameters are badly constrained (as ex- 
pected for a non significant excess) and discontinuity in the \ 2 
curves prevents us from evaluating the errors. Because it is a gen- 
eral property of QPOs that their RMS amplitude increases with en- 
ergy, we have also restricted the energy range from 4 to 20 keV, but 
failed to detect any significant excesses. We conclude that the 1330 
Hz QPO was likely a statist ical artifact, a possibility also implied 
by Ivan Straaten et all J2000j) . A summary plot of the analysis per- 
formed is presented in Fig(3] We also present the X-ray light curve 
corresponding to the ObsID, indicating no anomalies in the source 
behavior along the observation. 

2.2 Average properties of kHz QPOs 

Following on iBarret et"al] d2006f) . we can identify lower and up- 
per QPOs based on their position in the diagrams presented in Fig 
[2] One can then align QPOs of similar type (either a lower or an 
upper) using a shift-and-add technique dMendez et alJI 199*81) . This 
allows us to obtain a better description of their average properties; 
quality factor and RMS amplitude. We have therefore aligned all 
the ObsID averaged PDS, containing a QPO with R > 3 (either a 
lower or an upper), within a 50 Hz interval. As stated above, our 
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Figure 3. Top : Light curve of the ObsID 40030-01-04-00 in the energy 
band 4.6 keV - 97.3 keV. Second from the top : PDS in the energy ban d 
4.6 keV - 97.3 keV with a 5 Hz resolution, as in lvan Straaten et alTfeoOOl) . 
Third and fourth from top : PDS in the energy bands 2 keV - 20 keV, 4 keV 
- 20 keV with a 5 Hz resol ution. The blue curve rep resents the 1330 Hz 
QPO with parameters from Ivan Straaten et al.l l2OO0l) and the green curve 
re presents our best fit after i nitialization of the fit parameters with the values 
of lvan Straaten et alj i200ol) . 



conservative threshold ensures that we do add real QPOs. Through 
extensive simulations of PDS with QPO parameters appropriate for 
4U06 14+091 and statistics comparable to the real data, we have 
checked that with our procedure, which leads to averaging a large 
number of 1 second PDS (several thousands), the QPO parameters 
(quality factor and RMS amplitude) so recovered are not biased by 
any statistical fluctuations of the signal around the mean QPO pro- 
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Figure 4. Quality factor (left) and RMS amplitude (right) versus frequency after grouping the ObsIDs with QPOs of similar type and frequency (within a 50 
Hz interval). Red filled squares are for lower QPO. Blue filled circles are for upper QPO. An abrupt drop of coherence of the lower QPO around 700 Hz is 
now revealed. 



file. The results are presented in Fig. [4] This figure shows for the 
first time that the quality factor of the lower QPO starts by increas- 
ing with increasing frequency and then drops when it reaches a fre- 
quency around 700 Hz (the last point for the lower kHz is obtained 
after aligning all the PDS with an identified upper QPO above 1 100 
Hz). The maximum frequency averaged value of the quality fac- 
tor of the lower QPO is only about 25, while in one observation it 
reaches 60 albeit with large error bars. The low values reported can 
be explained in part by the fact that no correction for the frequency 
drift is applied within the ObsIDs. At the same time, the quality 
factor of the upper QPOs increases steadily. The behavior of the 
RMS amplitude of the lower QPO is consistent with other sources: 
it increases, saturates and then decreases sharply with increasing 
frequency. On the other hand, for the upper QPO, its RMS ampli- 
tude decreases up to a minimum around 1 100 Hz, after which a sec- 
ond maximum is observed at 1 150 Hz. Although the minimum was 
not co mpletely sampled by the data presented in lvan Straaten et all 
< |2002|) . a similar trend could be inferred. 



2.3 A constant frequency difference of twin QPOs 

As said above, in 24 ObsIDs, we detect simultaneous twin QPOs. 
Their frequency difference is plotted against the frequency of 
the lower QPO in Fig [5] This figure shows that albeit with 
large error bars, the frequency difference is consistent with be- 
ing cc4istani^_around320Hz. A similar conclusion was reached 
by Ivan Straaten et al.l d2000h using a reduced data set, but with a 
rather large scatter, with values ranging from 250 and 380 Hz. This 
value is indeed sig nificantly different from the neutron star spin fre- 
quency at 414 Hz dStrohmaver et alj|2008l) . Note also that the drop 
towards higher frequencies seen in other sources is not observed in 
4U 06 14+ 09, and our result is consis tent with findings of lFord et al.l 
dl997l) and lvan Straaten et alj feOOd) . 
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Figure 5. Frequency difference of the 24 simultaneous twin QPOs detected 
in our analysis. The frequency difference is consistent with being constant 
around 320 Hz. 



2.4 On the distribution of twin QPO frequencies 

We have started to investigate how the dependence of the QPO 
quality factor (or width) and RMS amplitude with frequency influ- 
ences their detectability in time intervals of fixed durations, hence 
affects the observed distribution of frequency ratios. This is an im- 
portant issue, because such distributions have been shown to be 
peaked around small integer ratios in some s ources, e.g. 3:2 in 
the case of Sco X -l dAbramowicz et al.ll2003h and 4U1636-536 
dTorok et alj|2008h . This has been used as an argu ment in favor 
of resonance based models (but see the discussion in lBelloni et al.l 
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Figure 6. Histograms of occurrence of single and twin upper QPOs show- 
ing a frequency gap between 1020 Hz and 1130 Hz in the distribution of 
frequencies. A frequency bin of 20 Hz has been used. The gap is indicated 
for twin upper QPOs as vertical dashed lines. 



d2005h ). Following this idea. lBarret & Boutelieil J2008I) have shown 
that in the case of 4U 1820-303, a gap in the frequency distribu- 
tion was present, together with a cluster of frequency ratio. They 
showed that the lack of twin QPOs within the gap could not be due 
to a lack of sensitivity for QPO detection, provided that the parame- 
ters of the QPOs (RMS and width) could be interpolated within the 
gap, using values measured before and after. Their result implied a 
sudden change of the QPO properties within the gap, most likely a 
loss of coherence of the upper QPO. 

As said above, there is a gap in the frequency distribution of 
the twin upper QPOs, between 1020 Hz and 1130 Hz. The his- 
tograms of occurrence of single and twin upper kHz QPOs are 
shown in Figure[6] 

Between 1020 Hz and 1 130 Hz, there is only one detection of a 
twin upper QPO and no single upper QPO detected. It is possible to 
evaluate the probability of having such a gap of 1 10 Hz width, with 
one QPO inside, assuming that the 24 twin upper QPO frequencies 
are uniformly distributed over their 250 Hz frequency span. This 
probability is less than 3.7 x 10~ 4 (3.4cr), giving us confidence 
that the gap is significant. As shown in Figure [6] single upper kHz 
QPOs are detected over a much wider range of frequencies and 
show a highly non uniform and peaked distribution of frequencies; 
single upper kHz QPOs are detected predominantly below the gap, 
but some are also detected above. 

There is no such a gap in the frequency distribution of the 
lower QPOs, 320 Hz below (there are more than 10 single lower 
QPOs). In fact, it is the frequency range in which both the quality 
factor and the amplitude of the lower QPO reaches a maximum; 
hence where lower QPOs are easy to detect. After recovering the 
upper QPO parameters within the gap through the shift-and-add 
technique, Figure [4] shows that the gap corresponds to the region 
where the RMS amplitude of the upper QPO reaches a minimum. 
This leads to the conclusion that the absence of twin QPOs within 
the gap is linked to a localized drop of the amplitude of the upper 
QPO, unlike in 4U 1820-303 for which the gap was more likely re- 



lated to a local minimum of the quality factor of the upper QPO. 
This result suggests some cautions when estimating frequency ra- 
tios based o n the detection of single lower OPOo Belloni et all 
120051 120071). Clearly in sensitivity limited observations (as is the 
case here, where the data are analyzed over comparable integration 
times), the histogram of ratios computed from simultaneous twin 
QPOs will differ from the one computed from the distribution of 
single lower QPOs. The study of the histogram of frequencies and 
frequency ratios, resulting from the frequency dependency of the 
QPO RMS and width (see Figure[4}, will be the subject of a forth- 
coming paper. 



3 CONCLUSIONS 

The main results of our systematic analysis of all archival RXTE 
data for 4U 0614+09 are: 

• We do not confirm the previous claim of a QPO at 1330 Hz. 
This is based on a thorough reanalysis of the observation from 
which the QPO was reported, and the fact that in our analysis, the 
highest frequency detected is at ~ 1220 Hz. This value is fully 
consistent with maximum frequencies observed in similar systems. 

• We observe for the first time a drop of the quality factor of the 
lower QPO. Such a drop has been interpreted as being related to the 
oscillating region, crossing the innermost stable circular orbit, and 
our detection is consistent with that idea. 

• The frequency difference between the lower and upper kHz 
QPOs is consistent with being constant, around 320 Hz, a value 
significantly different from the neutron star spin frequency (414 
Hz) or half its value. 

• If the drop in quality factor for the lower kHz QPO is due to 
the ISCO, one can estimate the orbital frequency there, given the 
constant frequency difference. To do this we compute the maxi- 
mum frequency of the lower QPO by extrapolating its quality fac- 
tor to zero. This yields a limiting frequency of ~ 920 — 930 Hz, 
corresponding to an orbital frequency of 1250 Hz at the innermost 
stable orbit (visco)- As a safety check, we note that the maxi- 
mum QPO frequency detected for the upper QPO is indeed lower 
than vi sco- From this we can estimate the mass of the NS fol- 
lowing the equation ^ « 2200 Hz x (1 + 0.75j)/visco where 
j = ~ 0-1 — 0.2 is the dimensionless angular momentum 
of the star. This leads to a gravitational mass for the neutron star 
of 1.9 Mq, i.e. a relatively massive NS, but still consistent with re- 
alistic modern equations of state, which predict maxi mum masses 
for slowly rotating ( j < 1) stars of ~ 1.8 - 2.3 M p, dAkmal et all 
1 19981 ; lLattimer & PrakastfeOOll ; iKlahn et alj|2006h . 

• 4U 0614+09 is the second object in which a gap in the fre- 
quency distribution of the upper QPO is observed, after 4U1820- 
303. This gap is likely associated with a local minimum of the RMS 
amplitude of the upper QPO. This is the first time such a minimum 
is so clearly identified, a lthough it was suggested in previous work 
( Ivan Straaten et alj2002l) . Similar trends will be searched in similar 
systems, because fr equency gaps are present in the data analyzed by 
iBarret et all j2006h in sources like 4U1636-536. This result moti- 
vates caution when estimating frequency ratios based on the detec- 
tion of single lower QPOs. 



2 and computing the frequency of the upper QPO through the linear func- 
tion that links the two QPO frequencies. 
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We note two additional consequences of our work. First, our con- 
clusion that the proposed 1330 Hz QPO in 4U 0614+09 is not likely 
to be real means that no kHz QPO source has a confirmed fre- 
quency significantly above ~ 1200 Hz. This is surprising a priori, 
because for a neutron star of canonical mass M = 1.4 M© and 
a spin parameter of j = 0.1, the orbital frequency at the ISCO is 
1690 Hz. Such frequencies should have been detected with RXTE, 
hence the wide gap relative to the maximum actually observed sug- 
gests a physical cause. One possibility is that these stars tend to 
have higher masses, consistent with accretion of several tenths of a 
solar mass, but other explanations should also be explored. 

Secon d, 4U 0614+09 is an important source to test the sug- 
gestion by iBarret et all d2005l . 120061 12007|) that the rapid drop in 
coherence of the lower QPO is caused by the approach of the oscil- 
lating region to the ISCO. In their proposal, the orbital frequency 
at the ISCO can be very roughly estimated by adding the separa- 
tion frequency Ai; between the twin kHz QPOs to the maximum 
frequency seen from the lower QPO after a sharp drop in qual- 
ity factor. If the 1330 Hz QPO were real, it would contradict this 
interpretation because the maximum inferred lower frequency is 
~ 900 Hz and the separation is ~ 320 Hz. However, the maximum 
significant upper QPO frequency is actually ~ 1220 Hz, which is 
consistent with the ISCO interpretation. This is encouraging, but 
other sources must also be analyzed carefully to look for potential 
disproofs of our hypothesis. 
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Table 1. The parameters of the single detected QPOs together with information concerning the observations: the ObsID name (ObsID), the date of the 
observation (Date), the start time of the observation (Time, UTC), the observation duration (T t, s ), the source count rate divided by the number of active PCU 
units (Rate), the background count rate (also divided by the number of active PCU units) (Bkg), the frequency (u), Full Width Half Maximum (FWHM), 
fractional RMS amplitude (RMS, %) and the R ratio, which is computed as the Lorentzian integrated power divided by its Icr error (the larger R, the more 
significant the QPO: R~ 3 corresponds to a ~ 6<r (single trial) excess in the power spectrum, see text for details). 



ObsID 




Date 


Time 


Tobs 


Rate 


Bk 


g 


V 


FWHM 


RMS 


R 


10072-02 


■02-00 


03- 


16 


1996 


09 


10 


3456 





92.8 


19 


.1 


725.0 ± 9.3 


108.3 ± 24.4 


16.3 ± 1.4 


6.0 


10072-02 


■02-00 


03- 


■16- 


■1996 


10 


45 


1120 





95.4 


19 


1 


753.3 ± 26.0 


154.9 ± 55.1 


17.5 ± 2.4 


3.6 


10095-01 


■01-00 


04- 


-22- 


■1996 


20 


56 


3440 


o 


89.0 


17 


.9 


565.4 ± 18.8 


133.1 ± 52.8 


15.9 ± 2.6 


3.1 


10095-01 


■01-00 


04- 


-22- 


1996 


22 


32 


3456 


o 


108.5 


17 


9 


601.0 ± 18.4 


128.6 ± 37.6 


14.1 ± 1.7 


1.2 


10095-01 


■02-00 


04- 


24- 


1996 




19 


2192 


o 


133.6 


20 


9 


582.8 ±2.1 


16.9 ± 7.3 


8.2 ± 1.1 


3.7 
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04- 


24- 
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18 


1 1 


1952 


Q 
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20 


9 
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109.7 ± 56.1 


11.9 ± 1.9 


3.1 
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04- 


94. 
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20 


58 
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o 
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20 


9 
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5.5 


10095-01 


■02-00 


04- 
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1996 


22 


34 
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20 


9 
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00 


10 
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01 


46 
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5.0 
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22 
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20 


9 
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08- 
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03 


30 
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o 
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114.0 ± 47.1 
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19 


18 
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17 


7 


587.2 ± 16.3 


144.5 ± 41.9 


17.5 ± 2.0 


1. 1 


20064-16 


■01-01 


08- 


15- 


1997 


00 


02 


3376 


o 


117.3 


19 


1 


636.0 ± 11.2 


97.3 ± 20.3 


12.7 ± 1.1 


5.7 


20074-01 


■01-00 


01- 


■10- 


1997 


10 


1 5 


3408 


o 


99.5 


17 


9 


585.5 ± 14.0 


134.3 ± 44.4 


15.3 ± 1.6 


4.8 


20074-01 


■02-00 


01- 


.25. 


1997 


1 3 


37 


3424 


o 


138.5 


21 


1 


641.7 ± 8.3 


109.6 ± 30.9 


13.1 ± 1.5 


4.3 


20074-01 


■03-00 


03- 


16 


1997 


22 


:2 3 


2960 


o 


124.4 


17 


9 


663.3 ± 11.8 


141.3 ± 35.8 


14.8 ± 1.5 


5.1 


20074-01 


■04-00 


05- 


■07- 


1997 


10 


12 


2240 


o 


158.9 


20 


.5 


739.3 ± 12.0 


109.1 ± 31.0 


11.8 ± 1.2 


1.9 


20074-01 


■05-00 


05- 


1 o- 


1997 


OS 


22 




n 
u 


138.6 


20 


.0 


655.1 ± 8.9 


110.7 ± 22.9 


13.0 ± 1.1 


6.0 


90074-01 


-06-00 


08- 


1 8- 


1997 


22 


28 


3440 





132.3 


20 


1 


574.3 ± 17.4 


151.1 ± 37.8 


12.6 ± 1.4 


4.6 


20074-01 


■07-00 


09- 


■03- 


1 QQ7 

lyy l 


2 1 


10 




n 

u 


126.6 


19 





636.4 ± 9.0 


117.5 ± 29.7 


14.2 ± 1.3 


5.4 


20074-01 


■08-00 


09- 


■06- 


1997 


[7 


;52 


2848 





126.6 


19 


,8 


621.2 ± 12.1 


134.2 ± 34.5 


14.3 ± 1.5 


4.7 


20074-01 


■10-00 


10- 


1 8- 


1997 


18 


19 


3456 


o 


109.1 


20 


1 


552.5 ± 19.2 


212.2 ± 56.9 


17.1 ± 2.0 


4.3 


30053-01 


■01-01 


10- 


■20- 


1998 


09 


04 


2496 


o 


83.9 


17 


6 


610.3 ± 18.8 


127.9 ± 46.1 


15.3 ± 2.2 


3.5 


30053-01 


■01-02 


10- 


20- 


1998 


13 


:57 


2256 


o 


85.0 


16 


9 


650.9 ± 13.9 


120.7 ± 46.5 


16.2 ± 2.2 


3.7 


30053-01 


■01-02 


10- 


■20- 


1998 




:28 


2976 


o 


86.7 


16 


9 


657.2 ± 16.8 


126.8 ± 42.0 


14.8 ± 1.9 


1.0 


30053-01 


■01-02 


10- 


■20- 


1998 


17 


04 


3408 


o 


89.8 


16 


9 


671.0 ± 11.2 


107.2 ± 32.9 


14.7 ± 1.6 


4.6 


30053-01 


■02-00 


10- 


16 


1998 


06 


04 


2800 


o 


133.6 


17 


9 


709.3 ± 3.3 


43.9 ± 6.8 


12.0 ± 0.7 


8.5 


30053-01 


■02-01 


10- 


J D- 


1 QQS 

1770 


02 


55 




n 

u 


128.9 


17 


.8 


653.6 ± 1.9 


23.5 ± 5.0 


9.6 ± 0.7 


6.5 


30053-01 


■02-02 


10- 


16 


1998 


18 


40 


3520 





133.9 


16 


.8 


682.7 ± 1.8 


28.3 ± 4.6 


10.3 ± 0.6 


8.6 


30053-01 


■02-05 


10- 




lyyo 


1 2 


1 7 


2448 


n 

u 


85.1 


17 


3 


653.7 ± 12.5 


88.7 ± 33.3 


17.6 ± 2.3 


3.8 


30054-01 


■01-00 


05- 


1 8- 


1998 


03 


58 


3424 





166.0 


21 


1 


781.8 ± 7.2 


82.4 ± 17.7 


10.4 ± 0.8 


6.3 


30054-01 


■01-00 


05- 


18- 


1998 


05 


34 


3456 





159.0 


21 


1 


736.7 ± 5.4 


90.4 ± 17.6 


12.7 ± 0.9 


6.8 


30054-01 


■01-02 


05- 


18- 


1998 


08 


:47 


3376 





157.5 


18 


3 


744.6 ± 9.6 


113.3 ± 24.1 


12.0 ± 0.9 


6.5 


30054-01 


■01-02 


05- 


18- 


■1998 


10 


:23 


3392 





155.4 


18 


3 


756.2 ± 7.5 


97.6 ± 21.1 


11.6 ± 0.9 


6.5 


30054-01 


■01-02 


05- 


18- 


1998 


11 


59 


1776 





155.7 


18 


3 


737.0 ±9.1 


112.5 ±22.6 


14.1 ± 1.1 


6.3 


30054-01 


■01-03 


05- 


19- 


1998 


03 


58 


3456 





151.4 


21 


6 


712.6 ± 9.4 


111.7 ±21.0 


12.5 ±0.9 


6.8 


30054-01 


■01-03 


05- 


19- 


1998 


05 


34 


1776 





154.9 


21 


6 


731.8 ± 10.6 


114.7 ±36.5 


13.3 ± 1.6 


1.3 


30054-01 


■01-04 


05- 


19- 


1998 


08 


46 


3440 





142.6 


18 


.2 


671.0 ±8.6 


120.7 ±23.6 


13.6 ± 1.1 


6,1 


30054-01 


■01-04 


05- 


19- 


■1998 


10 


:22 


2912 





139.7 


18 


.2 


647.1 ± 8.3 


124.7 ±27.6 


14.5 ± 1.3 


5.8 


30054-01 


■01-06 


05- 


17- 


■1998 


03 


58 


3440 





133.0 


21 


,1 


591.5 ± 9.8 


159.2 ±32.6 


16.5 ± 1.4 


6.0 


30054-01 


■01-06 


05- 


17- 


1998 


05 


34 


3456 





137.8 


21 


,1 


645.5 ± 11.4 


117.6 ±25.4 


12.5 ± 1.1 


5.9 


30054-01 


■02-00 


05- 


■22- 


1998 


02 


:23 


3184 





148.7 


20 


2 


727.8 ± 8.6 


114.8 ±22.7 


13.2 ± 1.0 


6.8 


30054-01 


■02-00 


05- 


22- 


1998 


03 


59 


3456 





151.6 


20 


.2 


738.8 ±8.3 


138.1 ±22.5 


14.8 ±0.9 


8.1 


30054-01 


■02-00 


05- 


■22- 


1998 


05 


35 


3456 





149.9 


20 


2 


737.2 ± 7.8 


100.0 ±21.0 


12.0 ±0.9 


6.1 


30054-01 


■02-01 


05- 


24- 


1998 


02 


:23 


3248 





134.1 


19 


9 


599.0 ±12.2 


116.6 ±27.9 


12.6 ± 1.2 


5.2 


30054-01 


■02-01 


05- 


24- 


1998 


03 


59 


3456 





124.1 


19 


9 


590.1 ± 12.5 


150.4 ±37.4 


14.6 ± 1.5 


4.8 


30054-01 


■02-02 


05- 


25- 


■1998 


00 


47 


2880 





147.3 


20 


,1 


694.5 ± 8.6 


111.1 ±22.6 


13.2 ± 1.0 


6,1 


30054-01 


■02-02 


05- 


-25- 


1998 


02 


:23 


3264 





149.9 


20 


,1 


735.6 ± 14.4 


162.1 ±30.8 


13.6 ± 1.1 


6.3 


30054-01 


■02-03 


05- 


26- 


1998 


02 


:23 


3328 





141.8 


19 


.5 


743.6 ± 8.8 


92.7 ± 22.3 


11.4 ± 1.0 


5.7 


30054-01 


■02-03 


05- 


26- 
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03 


59 
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19 
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11.0 ± 1.1 


1.8 
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Table 1. Continued 
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Table 2. The parameters of the twin QPOs together with information concerning the observations: the ObsID name (ObsID), the date of the observation (Date), 
the start time of the observation (Time, UTC), the observation duration (T t, s ), the source count rate divided by the number of active PCU units (Rate), the 
background count rate (also divided by the number of active PCU units) (Bkg), the frequency v ; u , Full Width Half Maximum (FWHM ; u ), fractional RMS 
amplitude (RMS ; u , %) and the R ; u ratio, which is computed as the Lorentzian integrated power divided by its lcr error, for the lower and upper QPO 
respectively. 
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10095-01-02-00 


04 


24 


1996 


14:57 


2080 





163 


2 


20 
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560.1 ±2.1 


13.3 ± 6.8 


6.2 ± 1.0 
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889.4 ± 10.5 


86.9 ± 34.2 


11.0 ± 1 
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8 


10095-01-02-00 


04 


24 


1996 


16:35 


1952 





166 
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21) 
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603.1 ± 6.2 


42.1 ± 12.7 


9.0 ± 1.1 
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932.4 ± 18.0 


84.0 ± 0.3 


9.4 ± 1 
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4 


10095-01-02-00 


04 


25 


1996 


04:58 


3456 





163 


5 


20 
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584.2 ± 4.0 


29.2 ± 10.8 


7.2 ± 0.9 
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80.3 ± 22.3 


10.3 ± 1 
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08 


06 
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10.3 ± 0.6 
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08 


07 


1996 


00:09 


3216 
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7 


18 


8 


628.2 ± 6.7 


54.9 ± 17.5 


9.4 ± 1.1 


4 


4 


944.4 ± 13.3 


108.0 ± 36.3 


11.3 ± 1 
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4 


4 


10095-01-03-000 


08 


07 


1996 


01:49 


2976 





146 


5 


18 


8 


582.9 ± 7.5 


36.8 ± 16.4 


7.1 ± 1.2 


3 





904.3 ± 12.2 


77.7 ± 26.4 


9.9 ± 1 


2 


4 





30053-01-02-00 


10 


16 


1998 


07:28 


3520 





133 


3 


17 


9 


709.8 ± 4.5 


60.9 ± 17.2 


11.8 ± 1.0 




8 


1017.0 ± 12.7 


107.2 ± 44.6 


10.1 ± 1 


4 


3 


7 


30053-01-02-00 


10 


16 


1998 


09:04 


3520 





120 


1 


17 


9 


648.9 ± 2.0 


27.1 ± 6.4 


9.3 ± 0.7 


6 


3 


977.4 ± 13.0 


118.4 ± 47.5 


11.1 ± 1 


.4 


4 





30053-01-02-00 


10 


16 


1998 


10:40 


3168 





110 


8 


17 


9 


580.5 ± 5.8 


36.7 ± 13.1 


8.4 ± 1.1 


3 


7 


893.3 ± 13.8 


126.4 ± 37.5 


13.2 ± 1 


4 


4 


7 


30053-01-02-02 


10 


16 


1998 


15:28 


2880 





117 


3 


16 


8 


572.3 ± 4.1 


34.4 ± 12.4 


8.8 ± 1.1 


4 


1 


896.6 ± 11.6 


72.1 ± 30.6 


9.2 ± 1 


4 


3 


3 


30053-01-02-04 


10 


16 


1998 


01:2(1 


2576 





126 


2 


19 


3 


640.0 ± 5.3 


43.2 ± 10.7 


9.9 ± 1.0 


5 


2 


954.8 ± 12.6 


59.5 ± 22.9 


8.8 ± 1 


3 


3 


5 


30056-01-01-02 


03 


07 


1998 


10:01 


3264 





1 11 


3 


17 


4 


741.5 ± 8.0 


84.9 ± 16.6 


11.3 ± 0.9 


6 


5 


1068.5 ± 11.8 


83.5 ± 30.4 


8.8 ± 1 


2 


3 


8 


30056-01-01-03 


0.; 


08 


1998 


10:20 


2144 





156 


2 


16 


9 


842.9 ± 19.8 


170.1 ± 73.0 


11.8 ± 1.7 
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5 


1144.4 ± 6.8 


67.9 ± 20.3 


10.2 ± 1 
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01 


1998 


15:25 


3168 
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1 


19 
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10.6 ± 0.8 
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311.1 ± 143.6 


13.7 ± 2 





3 


5 


40030-01-07-00 


12 
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1998 


22:02 


1520 





146 


2 


17 


8 


641.5 ± 1.9 


17.5 ± 5.3 


8.1 ± 0.9 


4 


7 


951.8 ± 19.4 


115.8 ± 48.5 


10.6 ± 1 


6 


3 


3 
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10 


21 


1999 


18:55 


3472 





136 


1 


19 


3 


800.4 ± 17.6 


109.6 ± 48.8 


10.2 ± 1.6 


3 


2 


1144.0 ± 9.4 


69.4 ± 25.4 


10.7 ± 1 


3 


4 


1 


80037-01-06-00 


03 


12 


2004 


05:28 


3104 







1 


30 





658.9 ± 2.9 


27.5 ± 6.4 


10.9 ± 1.0 


5 


7 


958.9 ± 11.4 


61.7 ± 25.5 


10.0 ± 1 


5 


3 


2 


80037-01-06-03 


03 


15 


2004 


07:30 


3248 





173 


3 


16 


5 


782.7 ± 12.8 


123.5 ± 41.4 


11.1 ± 1.2 


4 


5 


1130.6 ± 7.3 


45.8 ± 21.4 


7.3 ± 1 


1 


3 


2 


80037-01-09-01 


05 


04 


2004 


03:19 


3152 





162 


4 


16 


3 


591.5 ± 5.6 


35.8 ± 14.7 


8.2 ± 1.2 


3 


5 


921.3 ± 12.2 


75.4 ± 25.5 


10.0 ± 1 


3 


3 


9 


80051-02-02-00 


01 


24 


2004 


05:28 


2592 





136 


7 


16 


3 


696.7 ± 2.5 


33.9 ± 5.3 


11.0 ± 0.7 


8 


2 


1008.2 ± 12.9 


71.9 ± 29.0 


8.5 ± 1 


3 


3 


4 


80051-02-02-00 


01 


24 


2004 


07:02 


2800 





82. 


D 


16 


3 


673.4 ± 4.0 


30.1 ± 13.8 


10.3 ± 1.5 


3 


4 


992.2 ± 19.7 


75.5 ± 30.8 


10.8 ± 1 


.8 


3 





80051-02-02-00 


01 


24 


2004 


08:37 


2912 





78. 


7 


16 


3 


638.8 ± 8.7 


56.2 ± 20.6 


12.2 ± 1.6 


3 


8 


954.4 ± 12.7 


76.3 ± 34.9 


12.2 ± 1 


9 


3 


2 


90422-01-01-03 


10 


05 


2004 


18:07 


2816 





165 


9 


16 


4 


611.3 ± 9.3 


51.5 ± 14.3 


9.5 ± 1.1 


4 


2 


915.2 ± 22.9 


130.1 ± 44.7 


11.4 ± 1 


5 


3 


7 


92411-01-06-01 


10 


31 


2006 


19:13 


3200 





165 


1 


18 


7 


608.7 ± 3.5 


34.4 ± 11. 9 


9.7 ± 1.1 


4 


5 


928.8 ± 8.2 


69.7 ± 24.1 


10.8 ± 1 


3 


4 


2 



